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Abstract. Consider random ir-circulants A^,, with n — » oo 7 k = k(n) and whose input se- 
quence {ci/}/>o is independent with mean zero and variance one and sup„ n~ l E|«/| 2+l5 < 
oo for some S > 0. Under suitable restrictions on the sequence {k{ri)} n >i, we show that the 
limiting spectral distribution (LSD) of the empirical distribution of suitably scaled eigen- 
values exists and identify the limits. In particular, we prove the following: Suppose g > 1 
is fixed and pi is the smallest prime divisor of g. Suppose P g = I~n-i Ej where {Ej}\<j<g 
are i.i.d. exponential random variables with mean one. 

(i) If k s = -1 + sn where s = 1 if g — 1 and s = o(n p '^ i ) if g > 1, then the empirical 
spectral distribution of n^ l ^ 2 A/ l „ converges weakly in probability to UiP l J lg where Ui is 
uniformly distributed over the (2g)th roots of unity, independent of P g . 

(ii) If g > 2 and k g = 1 + sn with s = o(n pl_1 ) then the empirical spectral distribution 
of n- l/2 A k n converges weakly in probability to \J*>P g s where U2 is uniformly distributed 
over the unit circle in K 2 , independent of P g . 

On the other hand, if k > 2, k = with gcd(n,&) = 1, and the input is i.i.d. 

standard normal variables, then F n -\p_ At converges weakly in probability to the uniform 

distribution over the circle with center at (0, 0) and radius r = exp(E[log V^i]). 

We also show that when n = k 2 + 1 — > 00, and the input is i.i.d. with finite (2 + 5) 
moment, then the spectral radius, with appropriate scaling and centering, converges to the 
Gumbel distribution. 



For any (random) n x n matrix B, let H\(B), . . .,fi n (B) e C — R. 2 denote its eigenval- 
ues including multiplicities. Then the empirical spectral distribution (ESD) of B is the 
(random) distribution function on R 2 given by 



For a sequence of random « x n matrices {B n }„>\ if the corresponding ESDs Fb„ converge 
weakly (either almost surely or in probability) to a (nonrandom) distribution F in the space 
of probability measures on R 2 as n — » 00, then F is called the limiting spectral distribution 
(LSD) of {£„}„>!. See Bai (1999)[1], Bose and Sen (2007)[6] and Bose, Sen and Gan- 
gopadhyay (2009)[4] for description of several interesting situations where the LSD exists 
and can be explicitly specified. 

Another important quantity associated with a matrix is its spectral radius. For any matrix 
B, its spectral radius sp(B) is defined as 
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1. Introduction 



F B (x,y) = n l #\j : fij(B) € (-00, jc] x (-00, y], 1 < j < n). 




: /j. is an eigenvalue of B 
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where \z\ denotes the modulus of z e C. For classical random matrix models such as the 
Wigner matrix and i.i.d. matrix, the limiting distribution of an appropriately normalized 
spectral radius is known for the Gaussian entries (see, for example, Forrester(1993)[13], 
Johansson (2000)[19], Tracy and Widom (2000)[29] and, Johnstone (2001)[20]) which 
was later extended by Soshnikov[26, 27] to more general entries. 

Suppose a = {fl/};>o is a sequence of real numbers (called the input sequence). For 
positive integers k and n, define the nxn square matrix 



A k , n (g) = 



Qq U\ ... &n—\ 

Cl n -k @n-k+l • • • "n— &~1 
a n -2k a n -2k+\ ■■■ a n -2k-\ 



All subscripts appearing in the matrix entries above are calculated modulo n. Our conven- 
tion will be to start the row and column indices from zero. Thus, the Oth row of A^„(o) is 
{ao, fli, «2, a„_i).ForO< j < n-1, the (y'+l)-throw of Ak,„ is aright-circular shift of 
the y'-th row by k positions (equivalently, k mod n positions). We will write At „{a) = A^„ 
and it is said to be a k-circulant matrix. Note that Ai,„ is the well-known circulant ma- 
trix. Without loss of generality, k may always be reduced modulo n. Our goal is to study 
the LSD and the distributional limit of the spectral radius of suitably scaled fc-circulant 
matrices A^„(o) when the input sequence a = {fl/}/>o consists of i.i.d. random variables. 

1.1. Why study £-circulants? One of the usefulness of circulant matrix stems from its 
deep connection to Toeplitz matrix - while the former has an explicit and easy-to-state 
formula of its spectral decomposition, the spectral analysis of the latter is much harder 
and challenging in general. If the input {a/}/>o is square summable, then the circulant 
approximates the corresponding Toeplitz in various senses with the growing dimension. 
Indeed, this approximating property is exploited to obtain the LSD of the Toeplitz matrix 
as the dimension increases. See Gray (2006)[ i 4] for a recent and relatively easy account. 

When the input sequence is i.i.d. with positive variance, then it loses the square summa- 
bility. In that case, while the LSD of the (symmetric) circulant is normal (see Bose and 
Mitra (2002)[. ] and Massey, Miller and Sinsheimer (2007)[16]), the LSD of the (sym- 
metric) Toeplitz is nonnormal (see Bryc, Dembo and Jiang (2006) [7] and Hammond and 
Miller (2005)[15]) 

On the other hand, consider the random symmetric band Toeplitz matrix, where the 
banding parameter m, which essentially is a measure of the number of nonzero entries, 
satisfies m — > oo and m/n — > 0. Then again, its spectral distribution is approximated well 
by the corresponding banded symmetric circulant. See for example Kargin (2009) [2 1] and 
Bose and Basak (2009)[3]. Similarly, the LSD of the (n - l)-circulant was derived in Bose 
and Mitra (2002) [. ]) (who called it the reverse circulant matrix). This has been used in the 
study of symmetric band Hankel matrices. See Bose and Basak (2009)[3]. 

The circulant matrices are diagonalized by the Fourier matrix F = ((F Sj( )), F v = 
e lmst l n I V«, < s, t < n. Their eigenvalues are the discrete Fourier transform of the in- 
put sequence {a/}o</<« and are given by A, = 2"=o a/<?~ 2m '^", < t < n. The eigenvalues 
of the circulant matrices crop up crucially in time series analysis. For example, the peri- 
odogramof a sequence {a/}/>o is defined as n~'| Yj"~o aie 2m ^ n \ 2 , -L^J ^ V ^ L^J an d is 
a simple function of the eigenvalues of the corresponding circulant matrix. The study of 
the properties of periodogram is fundamental in the spectral analysis of time series. See 
for instance Fan and Yao (2003)[12]. The maximum of the perdiogram, in particular, has 
been studied in Mikosch (1999) [10]. 
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The Ai-circulant matrix and its block versions arise in many different areas of Mathe- 
matics and Statistics - from multi-level supersaturated design of experiment (Georgiou and 
Koukouvinos (2006) [ 1 7]) to spectra of De Bruijn graphs (Strok (1992)[28]) and (0, 1)- 
matrix solutions to A'" - J n (Wu, Jia and Li (2002) [31]) - just to name a few. See also the 
book by Davis (1979)[9] and the article by Pollock (2002)[23]. The fc-circulant matrices 
with random input sequence are examples of so called 'patterned' matrices. Deriving LSD 
for general patterned matrices has drawn significant attention in the recent literature. See 
for example the review article by Bai (1999)[ 1] or the more recent Bose and Sen (2008)[6] 
and also Bose, Sen and Gangopadhyay (2009)[4]). 

However, there does not seem to have been any studies of the general random A:-circulant 
either with respect to the LSD or with respect to the spectral radius. It seems natural to 
investigate these. The LSDs of the 1-circulant and 1-circulant with symmetry restriction, 
(n - 1) circulant are known. It seems interesting to investigate the possible LSDs that may 
arise from £-circulants. Likewise, the limit distributions of the spectral radius of circulant 
and the (« - l)-circulant are both Gumbel. It seems natural to ask what happens to the 
distributional limit of the spectral radius for general £-circulants. 

1.2. Main results and discussion. 

1.2.1. Limiting Spectral distributions. The LSDs for fc-circulant matrices are known for a 
few important special cases. If the input sequence {fl/};>o is i.i.d. with finite third moment, 
then the limit distribution of the circulant matrices (k = 1) is bivariate normal (Bose and 
Mitra (2002)[5]). For the symmetric circulant with i.i.d. input having finite second moment, 
the LSD is real normal, (Bose and Sen (2007)[6]). For the ^-circulant with k = n- 1, the 
LSD is the symmetric version of the positive square root of the exponential variable with 
mean one (Bose and Mitra (2002)[5]). 

Clearly, for many combinations of k and «, a lot of eigenvalues are zero. Later we 
provide a formula solution for the eigenvalues. From this, if k is prime and « = m x k 
where gcd(m, k) = 1, then is an eigenvalue with multiplicity (« - to). To avoid this 
degeneracy and to keep our exposition simple, we primarily restrict our attention to the 
case when gcd(£, n) - I. 

In general, the structure of the eigenvalues depend on the number theoretic relation 
between k and n and the LSD may vary widely. In particular, LSD is not 'continuous' in k. 
In fact, while the ESD of usual circulant matrices n~ l ^ 2 A\ A is bivariate normal, the ESD of 
2-circulant matrices n^^Aj^n for n large odd number looks like a solar ring (See Figure 1). 
The next theorem tells us that the radial component of the LSD of fc-circulants with k > 2 
is always degenerate, at least when the input sequence is i.i.d. normal, as long as k — n° m 
and gcd(£, «) = 1 . 

Theorem 1. Suppose {a?}/>o is an i.i.d. sequence of N(0, 1) random variables. Let k > 2 
be such that k — n o( - ' and n — > oo with gcd(n,k) — 1. Then F n -\n Akn converges weakly 
in probability to the uniform distribution over the circle with center at (0, 0) and radius 
r — exp(E[log V£]), E being an exponential random variable with mean one. 

Remark 1. Since — logE has the standard Gumbel distribution which has mean y where 
y ~ 0.57721 is the Euler-Mascheroni constant, it follows that r — e~ y l 2 =s 0.74930. 

In view of Theorem 1, it is natural to consider the case when k g - Q.{n) and gcd(£, n) = 1 
where g is a fixed integer. In the next two theorems, we consider two special cases of the 
above scenario, namely when n divides k g ± 1 . Consider the following assumption. 
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Figure 1. Eigenvalues of 100 realizations of n 1/,2 A^„ with a; ~ A^(0, 1) 
when (i) k = l,n = 901 (left) and (ii) k = 2,n = 901 (right). The color 
represents the height of the histogram - from red (high) to blue (low). 

Assumption I. The sequence {a/}/>o is independent with mean zero, variance one and for 
some 6 > 0, 

«-i 

supra" 1 V E\a,\ 2+S < oo. 
i=o 

We are now ready to state our main theorems on the existence of LSD. 

Theorem 2. Suppose {«/}/>o satisfies Assumption I. Fix g > 1 and let p\ be the smallest 
prime divisor of g. Suppose k g — -1 + sn where 5=1 if g — 1 and s — o(n Pl ~ l ) if 
g > 1. Then F n -\p. Akii converges weakly in probability to I/i(ITj-i ^/) * as n — > oo where 
{Ej}\<j< g are Ltd. exponentials with mean one and U\ is uniformly distributed over the 
(2g)th roots of unity, independent of '{E j)\<j< g . 

Theorem 3. Suppose {a/}/>o satisfies Assumption I. Fix g > 1 and let p\ be the smallest 
prime divisor of g. Suppose k 8 — l+sn where 5 = if g — 1 and s — o{n p ^ 1 ) if g > 1. Then 
F n -i/i Akii converges weakly in probability to U2(Y\ 8 j = i Ej)^ 2g as n — > oo where {Ej]\<j< g are 
i.i.d. exponentials with mean one and Ui is uniformly distributed over the unit circle in R 2 , 
independent of \E j}i<j< g . 

Remark 2. (1 ) Theorem 2 and Theorem 3 recover the LSDs of k-circulants for k — n — 1 
and k — 1 respectively. 

(2) While the radial coordinates of the LSD described in Theorem 2 and 3 are same, their 
angular coordinates differ. While one puts its mass only at discrete places e' 2n ^ 2s , 1 < j < 
2g on the unit circle, the other spreads its mass uniformly over the entire unit circle. See 
Figure 2. 

(3) The restriction on s — (k g ± l)/n in the above two theorems seems to be a natural one. 
Suppose g is a prime and so g — p\. In this case if s > n Pl , then k becomes greater than 
or equal to n violating the assumption that k < n. 

(4) We cannot expect similar LSDs to hold for more general cases like k 8 — +r + n, r > 1 
fixed. Compare Figure 2 and Figure 3. 
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Figure2. Eigenvalues of 20 realizations of n 1/,2 Aj„ witha; ~ Exp(l)— 1 
when (i) k = 1 1, k 3 = -1 + In (left) and (ii) k = 1 1, k 3 = 1 + 2n (right). 




Figure 3. Eigenvalues of 100 realizations of n ' 2 A^ n with a\ ~ N(Q, 1) 
when (i) k = 16, n = -3 + k 2 (left) and (ii) k = 16, n = 3 + k 2 (right). The 
color represents the height of the histogram - from red (high) to blue 
(low). 

1.2.2. Spectral radius. For the fc-circulant, first suppose that the input sequence is i.i.d. 
standard normal. When k = 1, it is easy to check that the modulus square of the eigenvalues 
are exponentials and they are independent of each other. Hence, the appropriately scaled 
and normalized spectral radius converges to the Gumbel distribution. But when the input 
sequence is i.i.d. but not necessary normal, that independence structure is lost. A careful 
use of Komlos-Major-Tusandi type sharp normal approximation results are needed to deal 
with this case. See Davis and Mikosch ( 1999)[ 10]. These approximations imply that the 
limit continues to be Gumbel. The spectral radius of the (n - 1 )-circulant is the same as that 
of the circulant and hence it has the same limit. See also Bryc and Sethuraman (2009)[8] 
who use the same approach for the symmetric circulant. 

Now let g — 2 and for further simplicity, assume that n — k 2 + 1 . If the input sequence 
is i.i.d. standard normal, then the modulus of the nonzero eigenvalues are independent and 
distributed according to {EiEo) 1 ^, where Ej,j= 1,2 are i.i.d. standard exponential. Thus, 
the behavior of the spectral radius is the same as that of the maxima of i.i.d variables each 
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distributed as {E\E2) l ^ A - This is governed by the tail behaviour of E\Ei. We deduce this 
tail behaviour via properties of Bessel functions and the limit again turns out to be Gumbel. 
Now, as suggested by the results of Davis and Mikosch (1999) [10], even when the input 
sequence is only assumed to be i.i.d. and not necessarily normal, with suitable moment 
condition, some kind of invariance principle holds and the same limit persists. We show 
that this is indeed the case. 

Theorem 4. Suppose {a/]/>o is an i.i.d. sequence of random variables with mean zero and 
variance 1 andE\ai\ y < oo for some y > 2. If n — k 2 + 1 then 

sp(n- l/2 A tn ) - d„ 



converges in distribution to the standard Gumbel as n — > oo where q — q(n) — |_f J and the 
normalizing constants c n and d„ can be taken as follows 

m (logn) 1 ^ 2 / lloglog«\ 1 7T 

(1) c„ = (81og»r 1/2 andd n = ( J. h + + rrTjlog-. 

Y2 \ 4 log« / 2(8 log ny' 1 2 

In the next section we state the basic eigenvalue formula for A:-circulant and develop 
some essential properties of the eigenvalues. In Section 3 and Section 4 we state and prove 
the results on LSD and the spectral radius respectively. An Appendix reproves the known 
eigenvalue formula for A:-circulant. 



2. Eigenvalues of the A:-circulant 

We first describe the eigenvalues of a &-circulant and prove some related auxiliary 
properties. The formula solution, in particular is already known, see for example Zhou 
(1996)[32]. We provide a more detailed analysis which we later use in our study of the 
LSD and the spectral radius. Let 

n-l 

(2) co = co„ := cos(2n/n) + /sin(27r/«), i 2 = -1 and A, = ^ a/aj' 1 , < t < n. 

1=0 

Remark 3. Note that {A,, < t < n) are eigenvalues of the usual circulant matrix A\^ n . 
Let p\ < p2 < ■ ■ ■ < p c be all the common prime factors of n and k. Then we may write, 

c c 

(3) n = n'Y[$ and k = kf\\p7- 

Here a q , f3 q > 1 and «', k' , p q are pairwise relatively prime. We will show that (« - «') 
eigenvalues of Ak, n are zero and «' eigenvalues are non-zero functions of a. 

To identify the non-zero eigenvalues of A^ n , we need some preparation. For any positive 
integer m, the set Z,„ has its usual meaning, that is, Z m — {0, 1, 2, ... , m - 1). We introduce 
the following family of sets 

(4) S (x) := {xk b mod ri : b > 0}, x € Z„. . 
We observe the following facts about the family of sets {S (x)} xe z n , ■ 

(I) Let g x = #S{x). We call g x the order of x. Note that go = 1. It is easy to see that 

S (x) = {xk h mod ri : < b < g x ). 
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An alternative description of g x , which we will use later extensively, is the following. For 

x G Z„<, let 

Ojc — {b > : b is an integer and xk b - x mod n'}. 
Then g x = min (9 v , that is, g A is the smallest positive integer b such that xk b = x mod «'. 

(II) The distinct sets from the collection {S(x)} xeZn , forms a partition of Z„>. To see this, 
first note that x e S (x) and hence U*ez < ^ 00 - Zb' ■ Now suppose S(i)nS (y) ^ 0. Then, 
jc^ 1 = jA:* 2 mod«' for some integers £>i,£>2 ^ 1- Multiplying both sides by k g *~ b ' we see 
that, x G S(y) so that, 5(x) c S(y). Hence, reversing the roles, S(x) = S(y). 

We call the distinct sets in {S(x)} xe z n , the eigenvalue partition of Z„- and denote the 
partitioning sets and their sizes by 

(5) P Q = {0}, P u ■ ■ ■ , Pi-\ and nj = #Pj, 0<j<{. 
Define 

(6) Uj ,:= Y\A*f«, i = 0, 1 €- 1. 

The following theorem provides the formula solution for the eigenvalues of A^ n . Since 
this is from a Chinese article which may not be easily accessible to all readers, we have 
provided a proof in the Appendix. 

Theorem 5 (Zhou (1996)[32]). The characteristic polynomial ofA^„ is given by 

i-\ 

(7) y L (A Kn ){A)=X n ~ n 'Y\(W-n^. 

j=o 

2.1. Some properties of the eigenvalue partition {Pj,0 < j < (}■ We collect some 
simple but useful properties about the eigenvalue partition in the following lemma. 

Lemma 1. (i) Let x,y e Z n >. Ifn' - to e S(y) for some to G S(x), then for every t e S(x), 
we have n' — t e S (y). 

(ii) Fix x 6 Z„'. Then g x divides g for every g G O x . Furthermore, g\ divides g x for each 
x 6 Z„-. 

( Hi) Suppose g divides g\. Set m :— gcd(k s — 1 , «')• Let X(g) and Y(g) be defined as 

(8) X(g) :— \x : x 6 Z„- and x has order gj, Y(g) :— \bn' ' /m : < b < mj. 
Then 

X(g) C Y(g), #Y(g) = m and \J X(h) = Y(g). 

h:h\g 

Proof, (i) Since t e S{x) - S (to), we can write t = tok 1 ' mod n' for some b > 0. Therefore, 
«'-/ = («'- f )A: fe mod n' € S («' - f ) = S (y). 

(ii) Fix ^ € O x . Since g A is the smallest element of O x , it follows that g x < g. Suppose, if 
possible, g = qg x + r where < r < g x . By the fact xg x = x mod n', it then follows that 

x = xk g mod n' — xk qSx+r mod n' = xk r mod «'. 

This implies that r e O x and r < g x which is a contradiction to the fact that g x is the 
smallest element in O x . Hence, we must have r — proving that g divides g\. 

Note that k g[ = 1 mod «', implying that xk 8 ' = x mod «'. Therefore g\ € (9 A proving 
the assertion. 
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(iii) Clearly, #Y(g) = m. Fix x e X(h) where h divides g. Then, xk g = xQj 1 ) 8 ' = x mod n', 
since g/h is a positive integer. Therefore «' divides x(k g - 1). So, n' jm divides x(k g - l)/m. 
But n' jm is relatively prime to (k 8 - l)jm and hence n' jm divides x. So, x = bn'/m for 
some integer b > 0. Since < x < «', we have < b < m, and x e T(g), proving 
Uh-.h\ g X(h) £ and in particular, X(g) c F(g). 

On the other hand, take < b < g. Then (bn'/m) k g = (bn'/m) mod«'. Hence, 
g e Obn'/m which implies, by part (ii) of the lemma, that g cn >/,„ divides g. Therefore, 
Y(g) c U/,.%, X(h) which completes the proof. □ □ 

Lemma 2. Let gi = <ff . . . qjj," where q\ < qz < ... < q m are primes. Define for 
1 < j < m, 

Lj := [qiAii 1 1 ' % : 1 < h < ■ ■ ■ < ij < m} 

and 

lj€Lj lj£Lj 

Then we have 

(i)#{x e Z„' : ^ < gi} = Gi - G 2 + G 3 - G 4 + • ••. 
f»JGi - G2 + G3-G4+ " <G,. 

Proof. Fix x 6 Z„<. By Lemma l(ii), g v divides g\ and hence we can write g x — q^' . . . qjj" 
where, < 77/, < jb for 1 < b < m. Since g x < g\, there is at least one b so that ijb < yi,. 
Suppose that exactly /2-many ^'s are equal to the corresponding y's where < h < m. To 
keep notation simple, we will assume that, - yt, 1 < b < h and ijb <Jb, h+1 < b < m. 
(i) Then x e Y(gijqb) for h + 1 < b < m and x i Y(g\/qi,) for 1 < b < h. So, x is counted 
(m - /z) times in G\. Similarly, x is counted y"^) times m &2, i^"^' 1 ) times in G 3 , and so on. 
Hence, total number of times x is counted in (Gj - G% + G3 — . . .) is 



/ m — h\ 




im — h\ 




fm — M 












\ 1 J 




\ 2 j 




i, 3 j 



(ii) Note that m-h > 1. Further, each element in the set {x e Z„> : g x < g\] is counted once 
in Gi - G% + G3 - . . . and (m - h) times in G\. The result follows immediately. □ 

2.2. Asymptotic negligibility of lower order elements. We will now consider the ele- 
ments in Z„, with order less than that of 1 e Z„, which has the highest order g\. We will 
need the proportion of such elements in Z„,. So, we define 

(9) v Kn ,:=—#{xe'Z n ,:gx<g\}. 

n' 

To derive the LSD in the special cases we have in mind, the asymptotic negligibility of 
Vk,n' turns out to be important. The following two lemmas establish upper bounds on 
and will be crucially used later. 

Lemma 3. (i) If gi = 2, then v k , n , = gcd(k - 1, «')/"'• 

(ii) Ifgi > 4 is even, and k g ^ 2 — -1 modn, then Vk, n < < 1 + ^ gcd(k g '^ b - l,n). 

b\g,,b>3 

(iii) If g\ > 2 and qi is the smallest prime divisor of g\, thenvt, n ' < 2n'^ l k gl ^ qi . 
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Proof. Part (i) is immediate from Lemma 2 which asserts that n'v^n 1 - #Y(l) = gcd(£ - 
l,n')- 

(ii) Fix x e Z„< with g x < g\. Since g x divides g\ and g x < g\, g x must be of the form g\jb 
for some integer b > 2 provided g\ jb is an integer. If b — 2, then xA: Sl ^ 2 = xk 8 * = x mod «' . 
But A;-'"/ 2 = -1 mod n' and so, xfc^ 2 = —x mod «'. Therefore, 2x — mod n' and x can be 
either or n' 12, provided, of course, n' 12 is an integer. But go — 1 < 2 < g\ /2 so x cannot 
be 0. So, there is at most one element in the set X(g\/2). Thus we have, 

#{x e Z„, :g x < gl } = #X(g l /2) + J] # f x 6 Z »' : ^ = Si lb) 

b\gu b>3 

= #X(gi/2)+ J] #X(gilb) 

%i,*>3 

<1+ #Y (8i/ b ) [by Lemma l(iii)] 

b\gu b>3 

= 1 + Y gcdO fl/fe - l,n') [by Lemma l(iii).] 

(iii) As in Lemma 2, let gi = q^qV ■ ■ ■ ql"' where q\ < qi < . . . < q m are primes. Then by 
Lemma 2, 

m 

n x v Kn , = Gi - G 2 + G 3 - G 4 + . . . < Gi = ^ gcd(F l/,ft - 1, «') 

2>=i 
m 

2>=i 

where the last inequality follows from the observation 

m m m m 

b=i b=\ b=\ b=\ 



Lemma 4. Lef fe anof c fee fwo fixed positive integers. Then for any integer k > 2, the 
following inequality holds in each of the four cases, 

gcd(Jfe*±l,Jk e ±l)<Jk gcd( *' e) + l. 

Proof. The assertion trivially follows if one of b and c divides other. So, we assume, 
without loss, that b < c and b does not divide c. Since, k c ± 1 = k c ~ b (k b + 1) + ± 1), 

we can write 

gcdOt* + l,k c + 1) = gcd(fc* + l,k c ~ b + 1). 

Similarly, 

gcd(yt* - l,k c + 1) = gcd(/fc fc - l,** -6 ± 1). 

Moreover, if we write c\ = c — \cjb\b, then by repeating the above step \cjb\ times, we 
can see that gcd(A:* ± l,k° ± 1) is equal to one of gcd(A:'' + 1, k c ' + 1). Now if c\ divides b, 
then gcd(fe, c) = c\ and we are done. Otherwise, we can now repeat the whole argument 
with b — c\ and c = b to deduce that gcd(A:'' + l,k c ' + 1) is one of gcdOt'' 1 + 1,1^ + 1) 
where b\=b — \ blc\\c\. We continue in the similar fashion by reducing each time one of 
the two exponents of k in the gcd and the lemma follows once we recall Euclid's recursive 
algorithm for computing the gcd of two numbers. □ 
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Lemma 5. (i)Fixg > 1. Supposed — — l+sn, n — > oo with s — 1 if g — 1 and s — o(n Pl ~ l ) 
if g > 1 where p\ is the smallest prime divisor of g. Then g\ — 2g for all but finitely many 
n and v^ n — > 0. 

(ii) Suppose k 8 = 1 + sn, g > I fixed, n — > oo w/r/j 5 = if g — 1 a«ii s — o(n Pi ~ y ) where p\ 
is the smallest prime divisor of g. Then g\ — g for all but finitely many n and v^ n — > 0. 

Proof, (i) First note that gcd(n,k) = 1 and therefore n' = n. When g — 1, it is easy to 
check that gi -2 and by Lemma 3(i), v^ ,, < 2/n. 

Now assume g > 1. Since Ic 2g = (sn - l) 2 = 1 mod n, g t divides 2g. Observe that 
gi + g — 2g/2 because k 8 — -1 mod «. 

If gi = 2g/b, where b divides g and b > 3, then by Lemma 4, 

gcd(F> - \,n) = gcd (k 2glh - l,(k 8 + \)l s) < gcd - l,k 8 + 1) < yt gcd < 2 ^ / '' *> + 1. 

Note that since gcd(2g/fo,g) divides g and gcd(2g /b,g) < 2g/b < g, wehave gcd(2g/b, g) < 
g/pi- Consequently, 

(10) gcdik 28 ^ - l,n) < fc s// " + 1 < (sn - l) Ul " + 1 = o(n), 

which is a contradiction to the fact that k 8 ' = 1 mod n which implies that gcd(&* -1, ri) = n. 
Hence, g\ = 2g. Now by Lemma 3(ii) it is enough to show that for any fixed b > 3 so that 
b divides g\, 

gcd(k 8 ^ b - l,n)/n = o(l) as n — > oo, 
which we have already proved in (10). 

(11) Again gcd(n, k) — 1 and «' = n. The case when g — 1 is trivial as then we have g x = 1 
for all x e Z„ and = 0. 

Since k 8 = 1 mod n, gi divides g. If gi < g, then g\ < glp\ which implies that 
k 8 ' < k 8 ^ Pl = (sn + l) 1 ^ 1 = o(n), which is a contradiction. Thus, g = g\. 
Now Lemma 3(iii) immediately yields, 

2k 8i/Pi 2(1 + sn) l/Pi 

Vk,n < < = 0(1). 

n n 

□ 

3. Proof of Theorem 1, Theorem 2 and Theorem 3 



3.1. Properties of eigenvalues of Gaussian circulant matrices. Suppose {a/}/>o are in- 
dependent, mean zero and variance one random variables. Fix n. For 1 < t < n, let us split 
A t into real and complex parts as A, — a hn + ib t n , that is, 




t < n. 
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For z e C, by z we mean, as usual, the complex conjugate of z- For all < t, f < n, the 
following identities can easily be verified using the above orthogonality relations 

B,(at,nh„) = 0, and E(a 2 n ) = E(fc 2 „) = „/2, 

A, = X n _ t , E(A,A,,) = nl(t + f = «), E(|^,| 2 ) = «. 
The following Lemma will be used in the proof of Theorem 2 and Theorem 3. 

Lemma 6. Fix k and n. Suppose that {fl;}o</< n are i.i.d. standard normal random variables. 
Recall the notations Pj and Tljfrom Section 2. Then 

(a) For every n, n _1 ' 2 a ti „, n~ x ' 2 b t ^, < t < n/2 are i.i.d. normal with mean zero and vari- 
ance 1 /2. Consequently, any subcollection {n ;i , Ylj, , . . .} of{Tlj}o<j<(, so that no member of 
the corresponding partition blocks [P 'j y ,P 'j 2 , ■ ■ ■} is a conjugate of any other, are mutually 
independent. 

(b) Suppose 1 < j <land <P j D (n - V j) = 0. Then all n-"' /2 Tlj are distributed as nj-fold 
product of i.i.d. random variables, each of which is distributed as E l ^ 2 U where E and U 
are independent random variables, E is exponential with mean one and U is uniform over 
the unit circle in R 2 . 

(c) Suppose 1 < j < I andPj = n - 9j and n/2 i <Pj. Then n-"J /2 Uj are distributed as 
(nj/2)-fold product of i.i.d. exponential random variables with mean one. 

Proof, (a) Being linear combinations of {fl/)o</< n , n _1 ^ 2 a t ,„,M _1/,2 fo () „, < t < n/2 are all 
jointly Gaussian. By (12), they have mean zero, variance 1 /2 and are independent. 

(b) By part (a) of the lemma, note that yr l l 2 A t = « _1 ^ 2 a, „ + in - x ' 2 b t „ is a complex normal 
random variable with mean zero and variance 1/2 for every < t < n and moreover, they 
are independent by the given restriction on Pj. The assertion follows by the observation 
that such a complex normal is same as E l ^ 2 U in distribution. 

(c) If t e Pj then n-te Pj too and t + n - t. Thus n~ l A t A„- t = n~ l {a 2 n + bj n ) which, by 
part (a), is distributed as Y/2 where Y is Chi-square with two degrees of freedom. Note that 
Y/2 has the same distribution as that of exponential random variable with mean one. The 
proof is complete once we observe that tij is necessarily even and the A t 's associated with 
Pj can be grouped into rij/2 disjoint pairs like above which are mutually independent. □ 

3.2. Proof of Theorem 1. Recall the notation Aj,t,Pj,nj and g x from Section 2. By 
Theorem 5, the eigenvalues of n~ l t 2 Ai c , n are given by 




where i = V^T. Fix any s > and < 6\ < 6 2 < 2n. Define 




Clearly, it is enough to prove that as n — > oo, 



(14) 
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Note that for a fixed positive integer C, we have 

c 

« _l ^ rij < rT x ^ #{1 < x < n : g x — u) 

\<j<C:nj<C u=2 
C 

< n~ l #{1 < x <n : xk" = x mod «} 

M=2 

C 



= n' 1 #{1 < x < n : x(k" - 1) = sn for some s > 1 ) 

M=2 

c 

< n _1 ^(Jk" - 1) < n l Ck c -> 0, as n -> oo. 



;/=2 



Therefore, if we define 



/-i 



^c= £ n Jt 

7=0: n/ <C 

then the above result combined with the fact that = {0} yields Nc/n — > 0. With C > 
(2k)/ (62 - #1), the left side of (14) can rewritten as 



1 r! f 27rs 

- V# \s: e [61,62], s = 



1,2,. ..,nj X 



(]~[|n- 1/2 ^|) 1/n; 'e(r-e,r + e) 



n — Nc 1 



7i,0 2 ],i = 1,...,«7 



7=0, n,>C 
£-1 



XI 



(02 -0l) 

2n 



— !— V nj xn] l #\s:- e (2ny l [9y 

(]^\n-^A t \) l/n} e(r-s,r + e)] + 0(^) 

rep, ) " 

0(C l )) X I ( f[ \n- l ' 2 A,\f h 6 (r - s, r + £ ) + O (^) 

( fl l^'^l) 1 '" e (r - s, r + e) + O^" 1 ) + O (^) 



7=0, nj>C 



(Ol-Ql) 

2n 



(15) 



(8 2 -0i) 1 
+ 



2tt n - yV c 



^ hjxI ([] |n _1/2 ^|) 1/nj <£(r-s,r + s) + 0(C l ) + O (-^ j 



7=0, «7>C V /eP, 



To show that the second term in the above expression converges to zero in L 1 , hence in 
probability, it remains to prove, 



(16) 



(]~[|«- 1/2 ^|) 1/ "^(r-£,r + e) 



te<P, 



is uniformly small for all j such that tij > C and for all but finitely many « if we take C 
sufficiently large. 

By Lemma 6, for each 1 < t < n, \n^ l ^ 2 A,\ 2 is an exponential random variable with mean 
one, and A, is independent of A,> if f + n — t and \A,\ = \Af\ otherwise. Let E, E\,Ei, ... be 
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i.i.d. exponential random variables with mean one. Observe that depending or whether or 
not T 3 ,- is conjugate to itself, (16) equals respectively, 

/ Hj/2 



( ]~[ E,) "' i(r-E,r + e) or P ([] VI,) t(r-s,r + s) 

, /=l ) V ,= 1 

The theorem now follows by letting first n — > oo and then C — > oo in (15) and by observing 
that Strong Law of Large Numbers implies that 



l c 



l/c 



| VI, — » r = exp(E[log VI]) almost surely, as C — » oo. 

v ,=i / 

□ 

3.3. Invariance Principle. For a set B c R rf ,rf > 1, let (SB)'' denote the '^-boundary' of 
the set B, that is, {dBf :- {y e R d : \\y - z\\ < rj for some z e dB). By <!>,/(■) we always 
mean the probability distribution of a cZ-dimensional standard normal vector. We drop 
the subscript 1 and write just O(-) to denote the distribution of a standard normal random 
variable. 

The proof of the following Lemma follows easily from Theorem 18.1, page 181 of 
Bhattacharya and Ranga Rao (1976)[2]. We omit the proof. 

Lemma 7. Let X\ , . . . , X m be W 1 -valued independent, mean zero random vectors and let 
V m — mT x Ylj=\ Cov(Xj) be positive-definite. Let G,„ be the distribution of m~ l ^ 2 T m (Xi + 
X2 + ■ ■ ■ + X m ), where T m is the symmetric, positive-definite matrix satisfying T^ t — V m l . If 
for some 5 > 0, E||Xy|| (2+<5) < 00 for each 1 < j < m, then there exist constants C, — Ci(d), 
i — 1,2 such that for any Borel set A C W 1 , 

m 

\G m (A) - <D rf (A)| < C im - s/2 (m- 1 V EHr^H^) + 2 sup O d ((dA)" - y), 

j=l veE rf 

< C,m- sl \A^ m {V m ))- {1 ^p 1+6 + 2 sup <5i(5A)" - y), 

yeW 

where p2+s = mT x YIJ=\ /t m in(^m) > is the smallest eigenvalue of V m , and 

r] = C 2 p2+sr^ 6/2 - 

3.4. Proof of Theorem 2. Since gcd(fc, ri) = 1, n' = n in Theorem 5 and hence there are 
no zero eigenvalues. By Lemma 5 (i), Vk^jn — > and hence the corresponding eigenvalues 
do not contribute to the LSD. It remains to consider only the eigenvalues corresponding to 
the sets Pj of size exactly equal to g\. From Lemma 5(i), g\ = 2g for n sufficiently large. 

n-l 

Recall the quantities rij = Wj, IT/ = n /e p,/l,, where A t = 2 < j < n. Also, 

for every integer t > 0, tk g — (-1 + sn)t = -t mod n, so that, A t and A„- t belong to same 
partition block S (t) = S(n — t), Thus each fly is a nonnegative real number. Let us define 

J„ = {0 < j < I : WPj = 2g], 

so that n = 2g#J n + nvi„. Since, u^„ — > 0, (#J„y l n — > 2g. Without any loss, we denote 
the index set of such j as J n = {1, 2, . . .#/„}. 

Let l,g, q 2 , . . .q 2s ~ 1 be all the (2g)fh roots of unity. Since rij = 2g for every j e J„, the 
eigenvalues corresponding to the set Pj are: 
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Hence, it suffices to consider only the empirical distribution of If'/" 5 as j varies over the 
index set /„: if this sequence of empirical distributions has a limiting distribution F, say, 
then the LSD of the original sequence rr x ^A^ n will be (r, 9) in polar coordinates where r is 
distributed according to F and 9 is distributed uniformly across all the 2g roots of unity and 
r and 9 are independent. With this in mind, and remembering the scaling -\/n, we consider 



F n ( X ) = {#j n y l J] i([n-«nj 5 < *) . 



Since the set of A values corresponding to any Pj is closed under conjugation, there exists 
a set J{j c Vj of size g such that 

Pj - \x : x e Mj or n - x e Jlj}. 

Combining each A, with its conjugate, and recalling the definition of and \b,^} in (1 1), 
we may write fly as 

First assume the random variables {fl/};>o are i.i.d. standard normal. Then by Lemma 6(c), 
F n is the usual empirical distribution of #J„ observations on (n*=i Ej) l/2g where {Ej}i<j< g 
are i.i.d. exponentials with mean one. Thus by Glivenko-Cantelli Lemma, this converges 
to the distribution of (Il* = i Ej) l ^ 2g . Though the variables involved in the empirical distri- 
bution form a triangular sequence, the convergence is still almost sure due to the specific 
bounded nature of the indicator functions involved. This may be proved easily by applying 
Hoeffding's inequality and Borel-Cantelli lemma. 

As mentioned earlier, all eigenvalues corresponding to any partition block Pj are all the 
(2g)th roots of the product n ; . Thus, the limit claimed in the statement of the theorem 
holds. So we have proved the result when the random variables {a/}/>o are i.i.d. standard 
normal. 

Now suppose that the variables {fl/}/>o are not necessarily normal. This case is tackled 
by normal approximation arguments similar to Bose and Mitra (2002)[5] who deal with 
the case k = n — 1 (and hence g — 1). We now sketch some of the main steps. 

The basic idea remains the same but in this general case, a technical complication arises 
as we need to control the Gaussian measure of the 77-boundaries of some non-convex sets 
once we apply the invariance lemma (Lemma 7). We overcome this difficulty by suitable 
compactness argument. 

We start by defining 



l/2g 

\ 7=1 



< X 



x e 



To show that the ESD converges to the required LSD in probability, we show that for every 
x > 0, 

E[F„(i)] -> F(x) and Var[F„(*)] -> 0. 

Note that for x > 0, 

#j„ 

E[F„(x)] = (#/„)-> J] Ffrrmj < x 2g ). 

./=! 
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Lemma 6 motivates using normal approximations. Towards using Lemma 7, define 2g 
dimensional random vectors 



X U = 2^ 



Cli cos 



2ntl 



a\ sin 



2ntl 
n 



teJlj\ < / < n, 1 < j < #/„. 



Note that 



E{X,j) = and n" 1 ^ Cov(X,j) = I 2g V /, 
i=\ 

Fix x > 0. Define the set A c R 2 # as 

g 

A :- {( v,v ; : 1 < y < g) : f] |2 ; «.v 2 + v 2 )| < x 2 *}. 
i=i 

Note that 

1=0 

We want to prove 

#j„ 

E[F n (x)] - 0> 2g (A) = (#/„)"' £ (P(n-*n^ < x 2 «) - 0> 2g (A)) ^ 0. 



/=! 



For that, it suffices to show that for every s > there exists N = N(s) such that for all 

n > N, 

n-l 

^n- 1,2 Y J X u eA)-^ 2g (A) 



sup 



< e. 



Fix e > 0. Find M] > large such that <t>([-Mi, M\\°) < £/(8g). By Assumption I, 
E(?T 1/2 a,,„) 2 = E(n- U2 b,, n ) 2 = 1/2 for any n > 1 and < t < n. Now by Chebyshev 
bound, we can find M 2 > such that for each « > 1 and for each < t < n, 

n\n- 1/2 a,J > M 2 ) < s/(8g) and P(|?T 1/2 Z?,,„| > M 2 ) < e/(Sg). 

Set M = max{Mi, M 2 ). Define the set B := {(xj,yj : 1 < j < g) e R 2 « : \xj\, \yj\ < M V j}. 
Then for all sufficiently large n, 

n-l n-l 
J (?T 1/2 Y € A) - 0> 2g (A) < P(«" 1/2 V € A n fi) - <D 2g (A n B) 

1=0 1=0 

We now apply Lemma 7 for A n B to obtain 

n-l 



+ s/2. 



K n ~ 1/2 X U eAnB )~ ®2g(A n B) 



1=0 



< 



C in - s/2 p 2+s + 2 sup <D 2g ((3(A n B)f - z) 

;eR 2 « 



where 



P2+S = P2+5 = sup n 1 V E||X /J || 2+ ' 5 and 77 = ?7(«) = C 2 p 2+(5 « 

0<l<n,jeJ„ ;=0 

Note that p 2+( s is uniformly bounded in n by Assumption I. 
It thus remains to show that 

sup 2g ((3(A n B))" - z) < e/8 



-5/2 
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for all sufficiently large n. Note that 

sup <b 2g ((d(A n B))' 1 -z)< sup I - zi) . . . 0(y g - Z2g)dx x . . . dy g 

< I rfxi . . . dy g . 

J(d(AnB)yi 

Finally note that d(A n B) is a compact (2g - l)-dimensional manifold which has zero 
measure under the 2g-dimensional Lebesgue measure. By compactness of d(A n B), we 
have 

(d(A n B)f i d(A r\B) as n -> 0, 

and the claim follows by Dominated Convergence Theorem. 

This proves that for x > 0, E[F„(x)] — » To show that Var[,F„(x)] — > 0, since the 

variables involved are all bounded, it is enough to show that 

n- 2 ^Cov(l(«^n 7 < x 2 «), l(n- g n r < x 2 «)) -> 0. 

M/' 

Along the lines of the proof used to show E[F n (x)] — > fX*), one may now extend the 
vectors with 2g coordinates defined above to ones with 4-g coordinates and proceed exactly 
as above to verify this. We omit the routine details. This completes the proof of Theorem 2. 
□ 

Remark 4. In view of Theorem 5, the above theorem can easily extended to yield an LSD 
has some positive mass at the origin. For example, fix g > 1 and a positive integer m. Also, 
fix m primes q\,qi, . . . ,q m and m positive integers /3\,B2, ■ ■ ■ ,B m . Suppose the sequences k 
and n tends to infinity such that 

(i) k — q\q2 ■ ■ -q m k and n — q^' q^~ . . .q^,"'h with k and h — > 00, 

(ii) k g — -1 + sh where s — o(n Pl ~ l ) — o(n Pl ■ ) where p\ is the smallest prime divisor 
ofg- 

Then F n -\n Akn converges weakly in probability to the distribution which has 1 - Tl'J =l qf J 
mass at zero, and the rest of the probability mass is distributed as f/i(II* = i Ej)^ 2s where 
U\ and {Ej}i<j< g are as in Theorem 2. 

3.5. Proof of Theorem 3. We will not present here the detailed proof of Theorem 3 but 
let us sketch the main idea. First of all, note that gcd(£, ri) = 1 under the given hypothesis. 
When g = 1, then k — 1 and the eigenvalue partition is the trivial partition which consists 
of only singletons and clearly the partition sets Pj, unlike the previous theorem, are not 
self-conjugate. 

For g > 2, by Lemma 5(ii), it follows that gi = g for « sufficiently large and Vk,„ — * 0. 
In this case also, the partition sets Vj are not necessarily self-conjugate. Indeed we will 
show that the number of indices j such that Pj is self-conjugate is asymptotically negligible 
compared to n. For that, we need to bound the cardinality of the following sets for 1 < b < 
81 = 8, 

W b := {0 < t < n : tk b = -t mod n] = {o < t < n : n\t{k b + 1)}. 

Note that to(b) := n/ gcd(rc, k b + 1) is the minimum element of Wi, and every other element 
of the set Wb is a multiple of to(b). Thus the cardinality of the set Wb can be bounded by 

#W h < n/t (b) = gcd(n,k b + 1). 
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Let us now estimate gcd(«, kr + 1). For 1 < b < g, 

gcd(n,k b + 1) < gcdt> ? - \,k b + 1) < k scd(gJ,) + 1 < K //J > + 1 = (1 + sn) 1/Pl + 1 = o(n), 
which implies 

n 1 J] #Wi = o(l) 

as desired. So, we can ignore the partition sets which are self-conjugate. 

Let J„ denote the set of all those indices j for which #Pj = g and Pj n (n - y 5 /) = 0. 
Without loss, we assume that J n = {1,2, . . .,#J„}. 

Let 1,£>, p 2 , . . . p* -1 be all the gth roots of unity. The eigenvalues corresponding to the 
set Pj, j € /„ are: 

n 1/g ,n 1/ V...n 1 V" 1 - 

7 7 J 

For j € J„, unlike the previous theorem n ; = Yltep ( a t,n + ibt.n) will t> e complex. 
Hence, we need to consider the empirical distribution: 

, #7„ g 

G„(x, y) = — 1 ( n J V" 1 < * + fy) , veR, 

" 7=1 r=\ 

where for two complex numbers w — x t + iy\ and z — x 2 + iyi, by w < z, we mean x\ < x 2 
and x 2 < y 2 . 

If {a/}/>o are i.i.d. ;V(0, 1), by Lemma 6, T^j g , j e !P ; are independent and each of them 

is distributed as ( FJjLi ^ as given in the statement of the theorem. This coupled 
with the fact that Q r ~ x U2 has the same distribution as that of U 2 for each 1 < r < g implies 
that {G„}„>i converges to the desired LSD (say G) as described in the theorem. 

When {fl/}/>o are not necessarily normals but only satisfy Assumption I, we show that 
EG n (x, y) — > G(x, y) and Var(G„(x, y)) — » using the same line of argument as given in the 
proof of Theorem 2. For that, we again define 2g-dimensional random vectors, 

Xij = 2 1/2 \ai cos |— j , a t sin ^— j : t e Pj j < / < n, 1 < j < #/,„ 

which satisfy 

n-l 

E(X,j) = and n _1 Y Cov(Xg) = / 2g V /, j. 
1=1 

Fix x,y € R. Define the set A c R 2 # as 



A : = Uxj, yj : 1 < j 

so that 



I \- '< v ; + /v ; ;| 



W=i 



< x + i v 



n-l 



(nJV- 1 < * + iy] = [n 112 ^ X,, e ^~'A\ 

1=0 

The rest of the proof can be completed following the proof of Theorem 2, once we realize 
that for each 1 < r < g, dig 8 r A) is again a (2g - l)-dimensional manifold which has 
zero measure under the 2g-dimensional Lebesgue measure. □ 
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4. Proof of Theorem 5 

We start by denning the gumbel distribution of parameter 6 > 0. 

Definition 1. A probability distribution is said to be Gumbel with parameter > if its 
cumulative distribution function is given by 

Ag(x) = exp{-flexp(-x)}, x e R. 

The special case when 6-1 is known as standard Gumbel distribution and its cumulative 
distribution function is simply denoted by AO- 
Lemma 8. Let E\ and E% be i.i.d. exponential random variables with mean one. Then 
(i) 

(17) K(x) := P(EiE 2 > x) = I exp(-y)exp(-xy _1 )c/y ~ tt 1/2 jc 1/4 exp(-2x 1/2 ) 

Jo 

as x — > 00. 

(ii) Let G be the distribution of (ExEi) ' . IfG, are i.i.d. random variables with the distri- 
bution G, and G™ := max {G, : 1 < t < n], then 

G (n) -d n d , 
> A. 

Cn 

where c n and d n are normalising constants which can be taken as follows 

i/o (logn) 1 ^ 2 / 1 loglog«\ 1 n 

(18) c n = (8\ogn)- l t 2 and d n = 1 + --p-g- + — -y^log-. 

Y2 \ 4 logn J /(slogn) 1 '- 2 

Proof, (i) Differentiating (17) twice, we get 
d 2 — f 00 

(19) T2 K W= ^ 2 exp(-y)exp(-xy- 1 )^, 
dx- Jo 

which implies that K satisfies the differential equation 

d 2 — — r°° t 

x—K(x) - K(x) = - (1 - xy- 2 ) exp(-(y + Jcy -1 ))^ 
dx l Jo 

(20) = exp(-(y + xy" )) =0, forx>0, 

with the boundary conditions /*T(0) = 1 and K(oa) = 0. 

From the theory of second order differential equations the only solution to (20) is 

K(x) = jrx 1/2 H\(2ix 1/2 ), x>0 

where i 2 = -1. The function H\ is given by (see Watson (1944)[30]) 

H\{x) = Jdx) + iY^x) 

where J\ and Y\ are order one Bessel functions of the first and second kind respectively. 

It also follows from the theory of the asymptotic properties of the Bessel functions J\ 
and Yi, that 

(21) K(x) x 7r 1/2 x 1/4 exp(-2x 1/2 ) as x 00. 
(ii) Now from (21), 

(22) G(x) = P{(E 1 E 2 ) 1/4 > x) x 7r 1/2 jcexp(-2x 2 ) as x 00. 
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By Proposition 1.1 and the development on pages 43 and 44 of Resnick (1996)[25], we 
need to show that, 

GO) = 0(x)(l - F#(x)) where lim 9(x) = 6 > 

X— >0O 

and, there exists some xq and a function / such that f(y) > fory > xq and such that / has 
an absolute continuous density with fix) — > as x — > oo so that 

(23) 1-F # 0) = exp(- f (1//Cy))dy), x>x . 
Moreover, a choice for the normalizing constants c„ and d„ is then given by 

(24) d* n = {\l{\-F # ))\n), <=/«). 
Then 

: * 

Towards this end, define for x > 1, 

(25) 0(x) x tt^V 2 , 1 - F # (x) = xexp(-2(x 2 - 1)), x> 1 = jc . 
To solve for /, taking log on both sides, 

(26) logx-2(x 2 -l) = - f J-rfy. 

J i /O0 

Taking derivative, 

1-2(2*) = --!- 
x /(x) 

or 

r(x) = — - x — as x — > oo. 

7 w 4x 2 - 1 4x 

Note that d* (to be obtained) will tend to oo as n — > oo. Hence 

c ;=/«)x(4<r 1 . 

We now proceed to obtain (the asymptotic form of) d*. Using the defining equation (24), 

(27) ^exp-^-'^n- 1 . 
Clearly, from the above, we may write 

/lognU/2 

where 5„ — > is a positive sequence to be appropriately chosen. Thus, again using (27), 
we obtain 

(log «)(<5 2 + 28 n ) - (i log log n + &) = 

where 

^ = 2--log2 + log(l + «. 

"Solving" the quadratic, and then using expansion yl+x — I + jx + 0(x 2 ) as x — > 0, we 
easily see that 



-2+ A /4 + 4(iloglog«+£,)/log« j 



^ loglog?7 +g„ ^ 
log n 



+ 



(log log n) 2 
(\ogn) 2 
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Hence 



' ± log log n + 
1 + 



+ 



(log log riy 
21og« j ~ \ (log«) 3 ^ 2 
Simplifying, and dropping appropriate small order terms, we see that 

G (n) -d„ d , 



where 



c„ = (81ogn) 1/2 and d„ 



(log/i) 1 ^ 2 / lloglogn 



V2 



1 + - 



4 log« / ' (8\ogn) l l 2 ^ 2 2 l ° gT> ' 



To convert the above convergence to standard Gumbel distribution, we use the following 
result of de Haan and Ferreira (2006) [1 1] [Theorem 1.1.2] which says that the following 
two statements are equivalent for any sequence of a„ > 0, b n of constants and any nonde- 
generate distribution function H. 

(i) For each continuity point x of H, 

lim G n (c n x + d„) = H(x), 

n— >°° 

(ii) For each x > continuity point of H~ (e~ 1 ' x ), 



lim 

t — >oo 



(1/(1 - G)Y (tx) - d[ t ] 



= H-\e~ xlx ) 



Now the relation A fl 1 (e ^) - A l (e l ' x ) = log 9 and a simple calculation yield that 

c n = c„, d„ = d n + c„ log^^V 2 ). 



4.1. Some preliminary lemmas. First of all, note that gcd(fc, ri) = 1 and hence n' = n. It 
is easy to check that g\ — A and 



be e : 



{0, n/2} if n is even 
{0} if n is odd. 



Thus the eigenvalue partition of {0, 1,2, ... ,n — 1} can be listed as 9*1,9*2, . . . ,f q , each of 
which is of size 4. Since each Pi, 1 < j < q is self-conjugate, we can find a set J[j c 9 > j of 
size 2 such that 

(28) Pj = {x : x e Jij or n - x e 
For any sequence of random variables b = {bi} t >o, define 

n-1 

(29) Pb,n{j)=n 2 X\ Yj b,t 



2ni 



, w = exp| — |, l<j<q. 



The next lemma helps us to go from bounded to unbounded entries. For each n > 1, define 



a triangular array of centered random variables {a] }o</<« by 

di - a" = fl//| fl ,|<„i/7 - Efl//| fl; |<„i/ r . 
Lemma 9 (Truncation). Assume E|a/| r < oo for some y > 2. Then, almost surely, 



max(/? fl ,„(;)) 1/4 - max06 5 ,„(;)) 1/4 = o(l). 

\<]<q \<j<q 
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Proof. Since Z"J oj" = for < t < n, it follows that B a , n (J) = Pa,n(J) where 

Cii-a [ = fl; + Ea;/| fl( |<„i/ r = fl;/| fl( |<„i/ r . 

By Borel-Cantelli lemma, with probability one, Yj'Zo \ a t\I\a,\>piy is finite and has only finitely 
many non-zero terms. Thus there exists an integer /V > 0, which may depend on the sample 
point, such that 

n-l n—i oo N 

(30) 2^ I 5 /"' ~~ a '\ = Zj l fl 'l / l«/l>» 1 "' ^ 2j ^ I M>' lh = X ^'IW' 1 ""- 



Consequently, if m > N, the left side of (30) is zero. Therefore, the terms of the two 
sequences {fl/}„,<;<„ and {a, } m <i<n are 
and the assertion follows immediately. 



sequences {«/},„<;<„ and {5, < " ) } m </<„ are identical almost surely for all sufficiently large n 



Lemma 10 (Bonferroni inequality). Let (£2, f, P) be a probability space and let B\, B2, . . . , B„ 

be events from T . Then for every integer m > 1, 

2m n 2m~l 

(3D J](-iy- l Sj,„<F(\jB,) < 2(-D ;_1 s *». 

7=1 J=l J=l 

l<il<l2<- ..<ij<r> l=\ 

Lemma 11. Fix x € R. Let E\,Ei, c„ and d n be as in Lemma 8. Let o~ 2 = n~ c , c > 0. Then 
there exists some positive constant K — K(x) such that 

P((£i£ 2 ) 1/4 > (1 + o- 2 „y 1/2 (c n x + d n )) <-, xel. 
v ' n 

Proof Since (1 +yT vl > 1 -y/lfoty > 0, 

P((£i£ 2 ) 1/4 > (1 + o- 2 „r U2 (c„x + d n )) < P((£,£ 2 ) 1/4 > (1 - o- 2 J2)(c n x + d n )) . 

Recall the representation 

P((£ , i£ 2 ) 1/4 > x) = 0(x)(l - F#(x)) as x -> 00. 

Note that (1 - cr 2 n l2){c n x + d„) = d* + (d„ - d*) + c„x - (c„x + d n )cr 2 l2 = d* + o x {Y) where 
we use the facts that c„ — > 0, (d„ - d*)/c„ = o(l) and d„ ~ ^J\ogn. The lemma now easily 
follows once we note that 1 - F#(d*) = re -1 , □ 

4.2. A strong invariance principle. We now state the normal approximation result and 
a suitable corollary that we need. For d > 1, and any distinct integers ■ ■ ■ , id, from 
{l,2,...,r^l}, define 

v 2d (l) = [cos |^ J , sin : 1 < j < rfj , I e Z„. 

Let ^x(-) denote the density of the 2c/-dimensional Gaussian vector having mean zero and 
covariance matrix S and let /2rf be the identity matrix of order 2c/. 
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Lemma 12 (Normal approximation, Davis and Mikosch ( 1999)[ 10]). Fix d > 1 andy > 2 
and let p n be the density function of 

n-\ 

2 1/ V 1/2 ^(a, + ^V / )v 2l/ (/), 

1=0 

where {Ni}i>o is a sequence of Ltd. N(0, 1) random variables, independent o/{a/};>o and 
cr 2 = Var(5o)s,r Ifn lc Inn < s\ < 1 with c — 1/2 - (1 - S)ly for arbitrarily small 5 > 0, 
then the relation 

p„(x) = tp (l+tT 2 )hd {x)(\ + s n ) with s„ -> 
holds uniformly for \\x\\ 3 = o d (n 1/2 - 1/r ), x e R M . 

Corollary 1. Let y > 2 and cr 2 — n~ c where c is as in Lemma 12. Let B c R M be a 
measurable set. Then 

I p n (x)dx- I ip (l+<T 2 )hi (x)dx <s„ I ip (1+a 2 )hd (x)dx + O d (exp(-n n )), 
Jb Jb Jb 

for some n > and uniformly over all the d-tuples of distinct integers 1 < ii < i% < . . . < 

i d < r^i. 

Proof. Set r = n a where < a < 1/2 - Using Lemma 12, we have, 

I p„(x)dx- I ip (l+CT i )hd (x)dx 
Jb Jb 

p„(x)dx- tp (l+cr 2 )l2ii {x)dx + J p„(x)dx+ ip (1+a .2 )hd (x)dx 

JBn{\M<r) JBn{\\x\\<r) J Bn{\\x\\>r) JBn{\\x\\>r) 

I V(\+aDi,S x )dx + I p„(x)dx+ <P(\+aj,)h d (x)dx = T] + T 2 + T 3 (say). 

Jfln{|MI<r) J{\\x\[>r] J{\\4\>r) 

Let v^J(0 denote the j-th coordinate of v 2 d(l), 1 < J < 2d. Then, using the normal tail 
bound, f(\N(0, cr 2 )| > x) < 2e-- x/cr for x > 0, 



T 2 



= I p n (x)dx = 

J{\M>r) 



n-i 



< 2d max I 

\<)<2d 



< 2d max I 

\<j<2d 



1=0 

n-\ ' 
2 l/2 w -l/2^ ( - / + anNl)v 0\^ > r/(2 d) 

1=0 

n x ' 2 ^ a/v^(/)| > r/(4 V2<f ) + 4<f exp ( - r«' /2 /(4 V2d)). 



Note that aivi (Z), < / < n are independent, have mean zero and variance at most one and 
are bounded by 2n^ r . Therefore, by applying Bernstein's inequality and simplifying, for 
some constant K > 0, 



|n _1/2 ^a/v^(/)| > r/4^/2d 



\ 1=0 



< exp(-Kr). 



Further, 



r 3 = I ^ (1+t7 2 )/ , rf (^)rfx < 4c/exp(-r/4fif). 

J{\\x\\>r) 



Combining the above estimates finishes the proof. 
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4.3. Proof of Theorem 4. First assume that n is even. Then k must be odd and S(n/2) = 
{n/2}. Thus with the previous notation, 

sp( t r ll2 A k ,„) = max(maxOS a ,„0-)) 1/4 ,kr 1/2 io|, |« 1/2 ^„ /2 |) . 

Since d q — > oo and c 9 — > 0, by Chebyshev inequality, we have 

sup P (\n~ ll2 A t \ > xc q + d q ) -> 0, for each x e R. 

0<f<« 

Thus finding the limiting distribution sp(«~ 1, ' 2 A / i j „) is asymptotically equivalent to finding 
the limiting distribution of nmx\<j< q (J3 a n {j)) 1 ^ . Clearly, this is also true if n is odd as that 
case is even simpler. 

Now, as in the proof of Theorem 2, first assume that {a;};>o are i.i.d. standard normal. 
Let {Ej)j>\ be i.i.d. standard exponentials. By Lemma 6, it easily follows that 

P( max(#,,„(f)) 1/4 > c q x + d q ) = ¥((E 2 j-iE 2j ) 1/4 > c q x + d q for some l<j< q). 

The Theorem then follows in this special case from Lemma 8. 

We now tackle the general case by using truncation of {fl/}/>o, Bonferroni's inequality 
and the strong normal approximation result given in the previous subsections. Fix x e R. 
For notational convenience, define 

Qf := P(max03 a+tr „;v,„O')) 1/4 > c q x + d q 

Qf := P(max(l + cr 2 n )(E 2M E 2j y /4 > c q x + d\ , 

where {Ni}i>o is a sequence of i.i.d. standard normals random variables. Our goal is to 
approximate Q " by the simpler quantity Qf. By Bonferroni's inequality, for all m > 1, 

2m 2m— 1 

(32) ^(_iy-is .„ < Q f) < ^ (-ly-'Sp,, 



7=1 7=1 

where 

S P' Yj P {^o-„N,n(h)) 114 > C q X + d q , . . . , (J3a +a - l)Njl {t > C q X + d q ) . 

l<h<t%<...<tj<q 

Similarly, we have 

2m 2m— I 

(33) ^(-ly-'r;,,, < of < £ (-ly-'r,;,,, 

7=1 7=1 

where 

Tu, = 

l<h<t2<...<tj<q 



2 P((l+cr2)(£ 2fl _ 1 £ 2 , l ) 1/4 > c g x+d q , (l+crlXE^E,,.) 1 ' 4 > c q x+d q ). 



Therefore, the difference between Qf and Qf can be bounded as follows: 



2m- 1 



(34) ^(-ly- 1 ^ - T ln ) - r 2m+1 ,„ < Qf - Qf < J] (-l) j -\S j,„ - TjJ + T 2mj „ 



7=1 7=1 
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for each m > 1. By independence and Lemma 11, there exists K = K(x) such that 

(35) Tjj , <[.)—<— foralln,;>l. 

Consequently, Hindoo lim sup n Tp, = 0. 

Now fix j> 1. Let us bound the difference between S p, and Tp,. Let defined in (28) 
be represented as 3\. t = {et, e' t )- For \- < h < tz < . . . < tj < q, define 



V2j(l) = 



1 2nle tl 



, sin 



2nle t 



, cos 



2nle' 



(2nle'\ 



. , cos 



(Inle'W 



Then, 



P((/?S+,r„iV,«(fl)) 1/4 > C q X +d g ,..., (fia+o- n Njtj)f IA > C q X + d q ) 

n-1 

= P(2 1/2 n- 1/2 Yffii + o-nNdv 2j (l) 6 B ( J } ), 



1=0 



where 
B 



U) 



{y e R 4 > : (y% +1 + yl +2 ) 1/4 (y 2 4l+3 + y 2 4t+4 ) l/4 > 2 1/2 (c 9 x + d q ),0 < t<j}. 



9 9 V 

By Corollary 1 and the fact +Ni — 2E\ , we deduce that uniformly over all the d- tuples 
1 < h < h < ■ ■ ■ < tj < q, 

n-\ 

(2 I/2 «- 1/2 Yjfli + <rJfi)v2jW e B ( J } ) - p((l + (r 2 ) I/2 (£ 2 ,,„_ 1 £ 2 ,,„) 1/4 > c q x + d q ,l<m< j) 



< e„P((l + ( r 2 ) 1/2 (£ 2/ „,_ 1 £ 2 ,„,) 1/4 > c q x + d q ,l< m < j) + C(exp(-«")). 
Therefore, as n — » oo, 

/«\ ^ 
(36) |S;,„ - r^l < e„r j; „ + I .jO(exp(-«")) < s„ — + o(l) -» 0, 

where 0{ ) and o(-) are uniform over j. Hence using (32), (33), (35) and (36), we have 

ill') ill') 

lim sup |<2j — Q\ I < lim sup r 2m+ i „ + lim sup r 2mi „ for each m > 1 . 
Letting ?7i — > oo, we conclude lim„ Qf - Qf = 0. Since by Lemma 8, 



max (£2,-1 £ 2j ) 1/4 = 0^((log«) 1/2 ) and at = n L ', 



it follows that 



and consequently, 



( 1 + crl) 1 ' 1 max^j^p^j) 1 ' 4 - d q d 



A 



maxi<j< q (B a+ITnNjl (j)y 14 - d q v 



A. 



In view of Lemma 9, it now suffices to show that 

max(/W,,iv,«0')) 1/4 - max (B^„(j)) 1/4 = o p (c q ). 

\<j<q \<j<q 

We use the basic inequality 

|kiZ2l - |wiw 2 || < (|zi| + |w 2 |)max{|zi - vvi|, \zi - w 2 |}, i\,zi,w\,W2 e C, 
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to obtain 



max(j3 s+a , iN ,„(j)) l/2 - max(J3- a ,„(j)) l/2 

l<j<q \<j<q 



< \M„{a + cr n N) + M n (a))M n (o- n N) 

< (2M„(a + cr„N) + M„(cr„N))M„(cr„N) 
where, for any sequence of random variables X = {X/}/>o, 



M n (X) :- max 

\<t<n 



n—1 



As a trivial consequence of Theorem 2.1 of Davis and Mikosch (1999)[10], we have 

M 2 n (o-„N) = O p (cr n logn) and M 2 (a + cr„N) = ; ,(log«). 
Together with cr„ = rT c ^ 2 they imply that 



From the inequality 



\<j<q \<j<q 



I y/yi - yjyi\ < . , 1. — p-rbi-^l, y\,yi>0 



it easily follows that 



maxOW^C/)) ' - max08 5 ,„(;)) 1/4 = o p («- c/5 ) = 0jJ (c 9 ). 

1<7<? \<j<q 

This completes the proof of Theorem 4. □ 

5. Concluding remarks and open problems 

To establish the LSD of £-circulants for more general subsequential choices of (k, n), 
a much more comprehensive study of the orbits of the translation operator acting on the 
ring Z„< by Tt(jc) = x£modn' is required. In particular, one may perhaps first establish 
an asymptotic negligibility criteria similar to that given in Lemma 5. Then, along the line 
similar to that of the proofs of Theorems 2 and 3 - first using the abundant independence 
structure among the eigenvalues of the £-circulant matrices when the input sequence is 
i.i.d. normals as given in Lemma 6 and then claiming universality through an appropriate 
use of the invariance principle. What particularly complicates matters is that in general 
there may be contributing classes of several sizes as opposed to only one (of size 2g or 
g) that we saw in Theorems 2 and 3 respectively. Thus it is also interesting to investigate 
whether we can select k = k{n) in a relatively simple way so that there exist finitely many 
positive integers h\,hi,..., h r , r > 1 with 

#{x e Z n , : g x = hj}/n' ^ Cj >0, l<j< r, 

where ci + . . . + c r = 1. In that case the LSD would be an attractive mixture distribution. 

Establishing the limit distribution of the spectral radius for general subsequential choices 
of (k, n) appears to be even more challenging. In fact, even under the set up of Theorems 2 
and 3, this seems to be a nontrivial problem. As a first step, one needs to find max-domain 
of attraction for (W S j = \ Ej)^ 2g where {£/}i<j< g are i.i.d. exponentials which requires a de- 
tailed understanding of the behaviour of P(TJ* =1 Ej > x) as x — > oo. When g > 2, we 
were unable to locate any results on this. Our preliminary investigation shows that this is 
fairly involved and we are currently working on this problem. Moreover, an extra layer of 
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difficulty arises while dealing with the spectral radius out of the fact that we can not im- 
mediately ignore some 'bad' classes of eigenvalues whose proportions are asymptotically 
negligible like we did while establishing the LSD. 
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Appendix 

Here we provide a proof of Theorem 5. Recall that for any two positive integers k and 
n, p\ < p2 < ■ ■ ■ < p c are all their common prime factors so that, 

c c 

n = n'Y\pf and k = k'\\p7 

where a q , f3 q >l and «', k', p q are pairwise relatively prime. Define 

(37) m :- max[/3 q /aq~\, [t] m j, := tk m mod b, b is a positive integer. 

\<q<c 

Let e m c i be a d x 1 vector whose only nonzero element is 1 at (m mod d)-th position, 
E m> d be the d x d matrix with ej m j, < j < d as its columns and for dummy symbols 
<!>o, Si,..., let A,„ £,d be a diagonal matrix as given below. 





(38) 

(39) 
(40) 

Note that 



,/xl 



E m ,d - [<?0,d e m j e2m,d ■ ■ ■ e(d-l)m,d] , 

A,„, M = diag[<S [0Lu) , S m ^ b , S [J]mJl 



Hd-l],„j, 



A ,m - diag [d mod b , 6 { mod h , 6 j mod h , 5 d _ x mod h ] . 
Lemma 13. Let n = (7r(0), ^(1), . . . , n(b - 1)) be a permutation of(0, \,...,b- 1). Let 

Ptt = [^K(0),b en(\),b ■ ■ ■ e„(b-\),b] ■ 

Then, P n is a permutation matrix and the (/, j)th element of P^E^bko^^Pn is given by 

(pTp i5, if(Uj) = (n- l (ktmodb),n- l (t)), 0<t<b 

(P n E k , b Aa Ab P„) u - | Q Qtherwise 

The proof is easy and we omit it. 

In what follows, x(A)(A) stands for the characteristic polynomial of the matrix A eval- 
uated at A but for ease of notation, we shall suppress the argument A and write simply 
X(A). 

Lemma 14. Let k and b be positive integers. Then 

(41) %(Ak,b) = x(E k ,bA),b,b) ■ 

where, 5 j — YX=q a i c ° jl ' ^ j <b, a> — cos(2n/b) + isin(2n/b), i 2 = —1. 
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Proof. Define the b x b permutation matrix 

Pb = 



4-i 

1 r 



Observe that for < j < b, the j-th row of A k , b can be written as a T pf where pf stands 
for jk-th power of Pb- From direct calculation, it is easy to verify that P/, = UDU* is a 
spectral decomposition of Pb where 

(42) D = diag(l, co, . . . , 

(43) U = [u ui ■■■ u h -i] with uj = b- x, \\,(j, co 2j , J b ~ l)j ), < j < b. 
Note that 6j = a T uj, < j < b. From easy computations, it now follows that 

so that, x (A k>b ) = x Ao,6,6), proving the lemma. □ 

Lemma 15. Let k and b be positive integers and, x — b/gcd(k, b). Let for dummy variables 
To, 71, 72, ...,76-1, 

F = (y , yi, y 2 , ••• , 7i>-i) ■ 

Then 

(44) X X T) = A b ~ x x(E Kx Xdiag(yQ modb , Jkmodb, Jlkmadb, ■ ■ ■ ,7(x-l)kmodb)) ■ 



Proof. Define the following matrices 

Bbxx = [eo.b eib eik,b 



e (x -x )Kb \ and P - [B B c ] 



where B c consists of those columns (in any order) of h that are not in B. This makes P a 
permutation matrix. 

Clearly, Ekj, — [B B ■ ■ ■ B] which is a b x b matrix of rank x, and we have 



t {E Kb T) = t (p T E Kb TP). 



Note that, 



P T E Kb TP = 



Q(b-x)xx Q(b-x)xx 



c 

Q(b-x)xb 



Q(b-x)Xx 



VP 



where, 



C 



C 

0(b-x)xb 



Vxix ■■■ /jr 



[B B c ] 



CB CB C 




= [/*/*■•■ Ix\ x diag(7o, 71, . . . , 7*-i). 

Clearly, the characteristic polynomial of P 1 ' Ekj,TP does not depend on CB C , explaining 
why we did not bother to specify the order of columns in B c . Thus we have, 

X(E k , b T) = x{P T E k j,rp) = A h - X % (CB). 
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It now remains to show that CB = E k , x x diag (y mod b, Jk mod *, Jik mod 7(x-i)k mod b)- 
Note that, the y'-th column of B is e_/jy,. So, y'-th column of CB is actually the {jk mod b)-t\\ 
column of C. Hence, {jk mod b)-\h column of C is jj k mo( j 6 ej k mo d t- So, 

CB = £jt^ X diag (yo mo d i» 7* mod fc, Jlk mod in • • • > TOc-l)* mod i>) 

and the Lemma is proved completely. □ 

Proof, of Theorem 5. We first prove the Theorem for A kn <. Since k and «' are relatively 
prime, by Lemma 14, 

X( A k,n>) = l{Ek,n'&Q,n',n>)- 

Get the sets So, S \ , ... to form a partition of {0, 1, . . . , n' — 1}, as in Section 2. 

Define the permutation n on the set Z„< by setting 7r(f) = 5,, < t < n' . This permutation 
7r automatically yields a permutation matrix P n as in Lemma 13. 

Consider the positions of 5 V for v e S / in the product P^E^AoxxPn- Let Nj-\ = 
Zin l^'l- We know, 5 y = {rjk x mod x > 0} for some integer rj. Thus, 

7T 1 (r^Jfe* -1 mod «') = Nj-i +t, 1 < f < «/ 

so that, position of S v for v = rjk'~ l mod 1 < t < rij in P^E kt ^Ao^P„ is given by 

Af/_i + 1 + 1, iV/_i + f) if, 1 < f < «/ 



iVj-i + 1, Nj-i+nj) 



if, f = « i 



(n~\rjk y mod n'^'Vi* 5 " 1 mod "')) : 
Hence, 

PlE Kn ,Ao, n ,, n ,P n = diag(Lo, Li, ...) 
where, for /' > 0, if n t ■ — 1 then Lj = \8 rj ] is a 1 x 1 matrix, and if rij > 1, then, 

... 5 ^,-1 > , 

rjk J mod «' 

<^mod„< ... 

S r . k mod n' . . . 



L> = 















0. 



r^A; J mod n' 

Clearly, x(Lf) = A n > - Ylj. Now the result follows from the identity 

x (£*,„' A ,„<,„') = ]~[ x(^) = n^" j - n ^ 

Now let us prove the results for the general case. Recall that n-n'y, TV =l p^''. Then, again 
using Lemma 14, 

X(Ak,n) = x(E ki „Ao M ). 
Recalling Equation 37, Lemma 14 and using Lemma 15 repeatedly, 

X(Ak,n) = X(E kM A ,n,n) 

= A»-"'x(E kj ,,A m+Pun ,) [forall;>0] 
= A"-"' x [E k>n , x diag ((5 [0]o „ , 6 Lvhji , <5 [2v]o „ , . . . , d [{n >-i )yhji )) [where y = n/n'] . 

Replacing Ao,„<,n' by diag (#[()]„„, 5[yj „ > £[2v] ,„> • • • > ^[(«'-i)v] ,„)' we can mimic the rest of the 
proof given for A k ^, to complete the proof in the general case. □ 
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